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ABSTRACT 

Seismograms  of  carefully  selected  large  earthquakes  recorded  on  the  high- 
gain,  long-period  network  were  used  to  visually  determine  the  amplitude  ratio 
of  20-second  Rayleigh  waves  in  the  ^ and  R3  phases.  Identification  of  R} 
was  aided  by  applying  a time-varying  processor  which  discriminates  heavily 
against  ground  motion  other  than  the  R3  signal.  Ratios  were  determined  in 
29  cases  representing  diveise  great  circles  over  the  earth;  and  the  average  Q 
estimated  for  the  20-second,  fundamental-mode  Rayleigh  waves  was  498.  This 
value  is  relevant  to  distance-amplitude  corrections  routinely  made  in  M 
calculations  and  is  significantly  greater  than  that  Q implied  by  the  current 
1.66‘logA  relation.  Our  data  did  not  reveal  any  significant  difference  in 
attenuation  over  oceanic  and  continental  structures  for  20-second  Rayleigh 
waves. 
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INTRODUCTION 


This  report  concerns  the  amplitude-distance  relation  of  Rayleigh  waves 
o'er  the  earth,  specifically  at  a period  of  20  seconds,  because  of  Its 
importance  in  the  routine  seismo logical  excercise  of  estimating  surface-wave 
magnitudes.  This  relation  can  be  approximated  analytically;  provided  that 
the  structure  is  laterally  invariant  and  that  frequencies  not  at  an  extremum 


of  the  group-velocity  curve  are  considered,  the  time-domain  amplitude  of 
seismic  surface  waves  of  a g!-/en  period  on  a sphere  is  given  by  the  relation 
(Sato,  1967): 

1 1 


or  taking  logarithms: 


sin  — 
r 


where  r and  rg  are  the  epicentral  distance  and  earth  radius  in  km  respectively, 
Q is  the  quality  factor  for  absorption,  U is  the  group  velocity  of  surface 
waves  with  period  T,  and  k is  a constant.  The  exact  rate  of  diminution  of 
amplitude  at  a given  period  is  mostly  affected  by  varying  Q values  since 
group  velocity  has  a small  range  of  values  for  a given  period  greater  than 
20  seconds,  but  it  is  important  in  assigning  surface-wave  magnitudes  to 
seismic  events  to  have  as  good  an  estimate  of  the  average  diminution  rate 
over  the  earth  as  possible.  Gutenberg  (1945)  found  that  observed  logarithms 
of  20-second,  Rayleigh-wave  amplitudes  between  15°  and  140°  distance  could  be 
fit  with  a straight  line  of  slope  1.656*logA,  where  A is  in  degrees;  that  is, 
a straight,  line  is  a good  approximation  to  the  terms  in  equation  (2)  over 
this  range.  Th  theoretical  relation  given  above  and  Gutenberg's  empirical 
approximation  to  it  are  shown  in  Figure  1 where  a range  of  Q values  has  been 
used  for  the  theoretical  calculations.  It  is  apparent  thaf  Gutenberg's 
relation  implies  a Q of  roughly  300  to  350  for  Rayleigh  waves  with  periods 
near  20  seconds.  Gutenberg's  early  relation  was  confirmed  by  the  Prague 
conference  (Vanek  et  al.,  1962)  which,  in  examining  many  proposed  amplitude- 
distance  relations,  found  1.66*logA  to  be  an  average  of  them;  since  then  the 
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1.66‘logA  relation  has  been  used  nearly  universally  in  computing  surface-wave 
magnitudes  of  seismic  events. 

In  some  recent  studies,  especially  with  underground  nuclear  explosion 
data,  it  has  been  found  that  Gutenberg's  rejation  was  improper  at  regional 
distances  (von  Seggern,  1970;  Evemden,  1971;  Basham,  1971)  and  that  a rela- 
tion using  roughly  1.0*logA  was  better.  These  problems  arose  simply  because 
Gutenberg's  formula  was  extrapolated  both  to  distances  and  to  periods  less 
than  it  was  intended  for  (Alewine,  1972).  However,  there  has  been  little 
questioning  of  its  correctness  beyond  regional  distances  except  for  a study 
(von  Seggern,  1975)  involving  a large  data  based  from  9 stations  in  the 
global  HGLP  (high-gain,  long-period)  network.  That  study  concludes  Rayleigh- 
wave  diminution  is  well  approximated  by  a 1.0*logA  relation  even  out  to  140°. 

This  relation  is  also  shown  on  Figure  1,  and  it  is  apparent  that  it  implies 
a Q on  the  order  of  1000  for  Rayleigh  waves  of  20  seconds. 

The  purpose  of  this  investigation  is  to  establish  a Q for  20-second 
Rayleigh  waves  which  will  reflect  a satisfactory  sampling  of  the  whole  earth 
and  therefore  be  appropriate  as  an  average  Q for  this  period.  Such  an  estimate 
wil.*.  be  relevant  to  improving  the  amplitude-distance  relation  for  Rayleigh-wave 
magnitude  determination  on  a routine  global  basis.  We  must  first  briefly 
review  the  estimates  already  published  on  this  parameter.  We  do  so  without 
discussing  the  possible  biases  in  the  studies  and  the  confidence  limits  on 
the  estimates. 

Arkhangelskaya  and  Fedorov  (1961)  found  Q ~ 340  for  22-24  sec  periods 
from  visual  analysis  of  R^  and  R£  phases  of  two  earthquakes  as  recorded  at 
the  Moscow  station.  However,  their  analysis  of  R£  through  R^  amplitudes  for 
the  same  events  gave  Q - 1000  for  these  periods.  Tryggvason  (1965)  used 
several  WWNSS  stations  and  one  explosion  on  Novaya  Zemlya  to  compute  spectral 
amplitudes  and  thus  Q for  a wide  range  of  frequencies;  his  result  for  20-sec 
Rayleigh  waves  was  Q = 1300.  Marshall  and  Carpenter  (1966)  estimated  Q ~ 400 
for  20-sec  Rayleigh  waves  using  visual  measurements  on  a suite  of  recordings 
in  the  Northern  Hemisphere  from  Novaya  Zemlya  explosions.  Burton  (1974) 
obtained  Q = 380  at  20  sec  from  Rayleigh-wave  spectral  calculations  at  many 
WWNSS  stations  for  Novaya  Zemlya  and  Lop  Nor  explosions.  Tsai  and  Aki  (1969) 
obtained  Q = 500  at  20  sec  from  the  same  type  of  analysis  on  an  earthquake. 


with  Q approaching  1000  for  periods  near  22  sec.  To  illustrate  the  elusive 
nature  of  the  parameter  we  seek,  20-sec  Q for  Rayleigh  waves,  we  point  to 
estimates  of  it  ranging  from  almost  100  to  2000  for  the  United  States 
from  spectral  calculations  in  Solomon  (1972),  Herrman  (1973),  and  Mitchell 
(1973). 

Certainly,  no  narrow  data  base  or  regional  study  such  as  those  mentioned 
here  will  provide  a Q estimate  -hat  will  be  representative  of  the  earth  as  a 
whole;  furthermore,  one  is  unsure  about  how  to  average  all  the  presently 
available  Q estimates.  Thus  a comprehensive  study  with  new  data  was  desirable. 


PROCEDURE 


Our  approach  to  the  problem  will  be  to  estimate  the  20-sec  Rayleigh-wave 
Q using  the  phases  R^  and  R^,  anti  we  will  do  this  with  visual  amplitude 
measurements.  The  appropriate  relation  is 


2ir2r 


Q = 


2ttt 


UT  (ln(A1/A3)  - j ln(l  + 


(3) 


This  is  derived  by  writing  equation  (1)  for  both  R^^  amplitude  A1  and  R3 

amplitude  A_,  taking  the  ratio,  and  solving  for  Q.  Note  that  the  geometric 

-1/2 

spreading  factor  sin  (x)  cancels  for  a complete  circuit  of  the  earth. 
Although  this  approach  is  the  same  as  Arkhangelskaya  and  Fedorov  used  in  a 
more  limited  study,  we  have  added  some  refinements  to  make  our  results  more 
robust.  The  validity  of  our  Q estimate  will  be  based  on:  1)  using  large 
events,  2)  examining  many  great  circle  paths  around  the  earth,  3)  isolating 
and  enhancing  the  R^  phases,  and  4)  certifying  that  no  signals  from  aftershocks 
or  other  events  are  masquerading  as  R^  phases.  Visual  measurement  is  deemed 
to  be  as  good  a procedure  as  spectral  measurements  in  this  case  because  only 
one  definite  period  is  being  considered,  Fourier  transformation  being  an 
unnecessary  step.  The  difference  between  a spectral  estimation  and  a visual 
estimation  is  merely  the  dispersion  factor  in  (1),  r_1^2,  which  comes  from 
the  stationary-phase  approximation  for  a travelling  surface  wave  (Sato,  1967). 
This  factor  disappears  in  the  case  of  spectrally-determined  amplitudes.  Also, 
in  the  case  of  spectral  estimation,  there  is  always  the  problem  of  window 
length;  for  the  20-second  amplitude  estimate  will  vary  with  this  length  in  an 
undefinable  manner.  In  conclusion  there  appears  to  be  no  real  advantage  to 
spectral  estimation  over  simple  visual  measurement. 

It  is  certain  here  that  we  are  mostly  concerned  with  amplitudes  given  by 

the  stationary-phase  approximation,  appropriate  to  periods  not  near  an  extremum 

of  the  group-velocity  curve,  rather  than  amplitudes  given  by  the  Airy-phase 

-1/3 

approximation  with  its  corresponding  r factor,  appropriate  to  periods  near 
an  extremum;  for  oceanic  group-velocity  curves  for  fundamental-mode  Rayleigh 
waves  have  no  extremum  near  20  seconds  and  only  a few  particular  continental 
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ones  do.  Paths  considered  in  this  study  are  50%  to  80%  oceanic,  and  so  the 
dispersion  factor  for 
great-circle  circuit. 


-1/2 

dispersion  factor  for  20-second  amplitudes  must  be  r for  most  of  each 
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DATA 


Data  for  this  study  was  taken  exclusively  from  stations  of  the  HGLP 
(high-gain,  long-period)  network.  The  locations  of  the  eleven  stations  jn 
this  network  are  given  in  Table  I.  Recordings  from  these  instruments  have 
magnifications  typicaily  ten  times  greater  than  WWNSS  long-period  instruments 
due  to  a shift  in  the  instrument  response  to  longer  periods,  rigid  environ- 
mental control,  and  installation  at  depth.  The  three-component  instrument 
outputs  are  also  digitally  recorded  at  a rate  of  one  sample  per  second. 

The  choice  of  events  was  governed  by  several  criteria.  Firstly,  they 
needed  to  be  large  enough  to  generate  visible  coda  to  4 hours  after  arrival 
of  R^;  this  gave  high  confidence  that  the  phase  R^  was  detectable  and  identi- 
fiable at  stations.  Secondly,  they  needed  to  be  small  enough  that  quality 
control,  timing,  and  calibration  of  the  digital  traces  could  be  carried  out 
by  use  of  the  paper  recordings.  These  criteria  translated  to  events  of 

t 6 or  Ms  2 6,  for  which  most  body-wave  phases  were  readable  but  most  R^ 
amplitudes  exceeded  the  available  width  on  the  paper  seismograms.  Thirdly, 
events  were  chosen  only  if  no  aftershocks  in  the  3-4  hour  period  after  R^ 
were  reported  on  available  epicenter  lists.  Fourthly,  a global  distribution 
of  epicenters  was  desired;  that  is,  m~ny  redundant  observations  due  to  nearby 
epicenters  was  to  be  avoided. 

Point  (3)  above  is  a crucial  requirement;  and  clearly  the  epicenter  lists 
must  be  complete  down  to  magnitude  4 if  we  are  to  be  sure  that  no  aftershocks 
are  responsible  for  erroneous  R^  picks  when  a R^/R^  ratio  of  roughly  100  is 
assumed.  However,  the  C-list  threshold  is  near  m^  z 5 globally  while  the 
LASA  and  NORSAR  thresholds  a'-e  m^  z 4 only  for  small  parts  of  the  globe.  We 
therefore  sought  to  eliminate  the  aftershock  uncertainty  by  checking  short- 
period  recordings  of  available  stations  closest  to  the  epicenters  chosen. 

For  each  event  we  verified  that  no  aftershocks  were  present  at  a time  which 
might  allow  an  R^  from  the  smaller  aftershock  to  be  mistaken  for  R^  of  the 
main  shock.  We  augmented  this  short-period  analysis  with  long-period  match 
filtering  of  time  windows  for  R^  using  the  R^  of  the  main  shock.  This  cross- 
correlation procedure  (Capon  et  al.,  1969)  should  definitely  identify  within 


TABLE  I 

Locations  of  Stations  in  the 
High-Gain,  Long-Period  Network 
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the  R^  window  any  R^  from  an  aftershock  or  from  an  unrelated  event  within  a 
few  hundreds  of  kilometers  from  the  epicenter  of  the  main  shock. 

fwelve  events  meeting  the  above  criteria  were  selected  for  processing 
and  are  listed  in  Table  II.  The  choice  of  stations  for  each  event  w&s  made 
on  the  basis  of  availability  and  data  quality;  over  one-half  of  the  possible 
vent— station  combinations  were  unable  to  be  processed  due  to  station  down 
time  on  one  or  more  components  or  inability  to  retrieve  the  digital  data  from 
the  tapes.  Some  other  combinations  were  eliminated  when  visual  analysis 
showed  that  coda  amplitude  to  be  below  normal  background  noise  level  within 
3 to  4 hours  after  R^;  these  noise  levels  varied  by  nearly  an  order  of  magni- 
tude. R^  itself  was  often  lower  in  these  cases  due  most  likely  to  radiation 
pattern.  So  for  these  cases  a combination  of  source  effects  and  station 
noise,  not  necessarLly  high  attenuation  for  R^,  made  them  unsuitable  for  this 
study. 
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PROCESSING 


The  three  components  of  the  36  qualifying  seismograms  were  calibrated 

relative  to  one  another  and  rotated  into  components  parallel  and  transverse 

to  the  great  circle  path,  with  the  vertical  component  unchanged.  Then  a 

Rayleigh-wave  enhancer  as  described  by  Simons  (1968)  was  applied  to  the  entire 

time  interval  from  the  start  of  R^  past  the  expected  time  of  arrival  of 

20-second  R^.  This  time-varying  signal  enhancer  operates  in  the  frequency 

domain  to  eliminate  harmonic  components  not  of  the  desired  Rayleigh-wave  type. 

It  discriminates  against  ground  motion  which  is  Rayleigh  waves  but  off  the 

4 

great  circle  path  by  a factor  equal  to  cos  a where  a is  the  difference 
between  the  observed  and  computed  back  azimuth.  Cases  where  R^  arrives  at 
angles  significantly  off  the  back  azimuth  are  not  unusual;  and  it  is  expected 
that  for  R^,  which  travels  an  additional  full  circuit  of  the  earth,  these 
cases  would  be  more  common  and  certainly  more  severe.  So,  even  though  we 
are  discriminating  against  signals  from  other  azimuths,  we  may  be  attenuating 
the  R^  phase  itself  somewhat;  and  so  our  Q estimates  will  be  biased  to  the 
low  side.  Also,  the  processor  discriminates  against  harmonic  components  not 

g 

having  the  Rayleigh-wave  particle  motion  by  a factor  equal  to  cos  8 where  6 
is  the  difference  between  the  observed  phase  angle  of  the  vertical  relative 
to  the  radial  recording  and  the  required  phase  angle  of  ir/2.  Since  the 
presence  of  considerable  noise  in  the  R^  window  should  cause  deviations  in 
this  angle,  we  expect  some  further  attenuation  of  the  R^  phase;  and  this  will 
further  bias  our  Q estimates  to  the  low  side. 


j 
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ANALYSIS 


Results  of  individual  analysis  of  the  36  processed  seismograms  are 
presented  in  Table  III.  Loth  and  R^  amplitudes  were  measured  on  t^e 
processor  outputs  for  the  vertical  component.  The  critical  sta^e  of  the 
analysis  was  the  identification  of  the  R3  phase  having  periods  near  20  seconds. 
This  identification  was  aided  by  checking  the  group  velocity  indicated  by  the 
20-second  R^-R^  time  delay  against  expected  group  velocity  based  on  the 
structure  of  the  complete  great  circle  path.  These  predicted  paths  are  shown 
in  Figure  2.  Percentages  of  oceanic  and  continental  structure  were  measured 
in  each  case;  and  using  3.8  km/sec  for  the  former  and  3.0  km/sec  for  the 
latter,  we  calculated  the  expected  travel  time  and  average  group  velocity  for 
each  case.  The  observed  group  velocities  are  listed  in  Table  III;  recognizing 
that  significant  departures  from  the  great  circle  path  might  occur  for  R^,  we 
accepted  as  valid  some  R^  arrivals  which  were  considerably  earlier  or  later 
than  expected,  where  observed  group  velocity  was  up  to  0.2  km/sec  different 
from  that  expected.  Identification  of  R^  arrivals  in  the  20-second  period 
range  was  also  frequently  aided  by  visibility  of  the  entire  R^  phase,  smoothly 
dispersed  from  a ionger-period  Airy  phase  down  through  shorter  periods  below 
20  seconds. 

The  20-second  R^  phases  were  graded  as  indicated  in  Table  III  from  0 to 
3 according  to  the  quality  or  confidence  of  the  observation.  Quality  "3" 
indicated  a 20-second  R^  which  is  unmistakable  at  even  a casual  glance;  if  so 
graded,  the  phase  is  most  often  identifiable  on  the  unprocessed  seismogram 
and  is  part  of  a plainly  dispersed  R3  wave  train.  Quality  "2"  indicates  an 
R^  which  is  identifiable  only  after  processing;  some  examination  of  the 
seismogram  is  required  to  establish  its  validity,  but  it  is  certainly  felt 
to  be  a correct  pick.  Quality  "1"  indicates  an  R^  of  uncertain  validity; 
these  phases  were  identified  only  after  careful  consideration  of  all  relevant 
data,  anl  there  is  no  clearly  dispersed  wave  tra'.n.  Quality  "0"  indicates 
an  R.  could  not  be  identified;  either  noise  obscured  it  or  its  angle  of 
approach  was  so  different  from  the  back  azimuth  to  the  epicenter  that  the 
time-varying  adaptive  processor  attenuated  it  severely.  An  example  of  each 
of  the  four  qualities  is  shown  in  Figure  3. 
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Great-circle  paths  of  36  signals  processed  for  Q determination 
A)  Events  1-7  , B)  Events  8-12. 


Figure  2 


RESULTS 


The  Q values  in  Table  III  were  calculated  according  to  equation  (3). 
Assuming  the  amplitude  ratios  A^/A^  to  have  log-normal  distribution,  the 
individual  1/Q  estimates  would  be  normally  distributed.  On  this  basis  the 
mean  and  90%  confidence  limits  for  1/Q  were  then  calculated  and  inverted, 
giving  Q = 441  and  408  < Q < 478.  These  confidence  intervals  are  not  meaning- 
ful except  to  indicate  the  scatter  of  the  data  because  we  have  already  explained 
that  our  average  Q estimate  will  be  biased  low  due  to  the  processor  attenuating 
Rj  more  than  R^.  In  order  to  establish  quantitatively  how  great  this  bias 
might  be,  we  reexamined  the  phases  on  those  recordings  having  quality  3 
identifications  of  R^.  In  10  of  the  12  cases,  R^  was  reliably  measurable  on 
the  unprocessed  trace  by  simply  correlating  the  motion  before  and  after 
processing;  in  other  words,  noise  was  little  or  no  hindrance.  With  R^  also 
remeasured  on  the  unprocessed  trace,  we  were  able  to  form  A^/A^  ratios  for 
the  raw  data;  these  ratios  were  uniformly  lower  than  these  presented  in 
Table  III  for  the  processed  data  by  a factor  of  between  one-fourth  and  one. 

The  average  factor  over  the  10  cases  was  .63;  assuming  this  is  a rough  empirical 
correction  factor  for  all  A^A.^  ratios  in  Table  III,  a corrected  set  of  29  ratios 
was  used  to  obtain  Q = 498.  We  regard  this  as  the  more  accurate  estimate. 

Since  we  have  graded  the  R^  identifications,  it  is  tempting  to  delete 
the  lesser  quality  R-^/R^  ratios,  flagged  as  1 in  Table  III,  and  recompute 
the  average  Q.  But,  noting  that  these  ratios  in  Table  III  result  in  the 
lowest  Q's  generally,  for  high  attenuation  paths  we  regard  these  as  valid 
estimates  which  should  be  equally  weighted  in  computing  the  mean  Q over  the 
globe.  For  those  seven  cases  where  no  R^  was  identified,  and  thus  where 
attenuation  might  be  considered  to  be  greatest,  an  upper  bound  for  R^  amplitude 
can  be  set  by  measuring  the  largest  20-second  cycle;  Q's  calculated  on  this 
basis  though  were  not  generally  low.  We  feel  that  in  these  seven  cases,  the 
inability  to  identify  R^  is  due  as  much  to  high  noise  or  to  propagation  far 
off  the  back  azimuth  as  to  low  Q;  and  we  refrain  from  using  these  noise 
amplitudes  to  make  a Q estimate. 

Since  we  have  already  estimated  the  percentage  of  oceanic  and  continental 
structure  in  each  path,  it  is  desirable  to  test  whether  our  data  indicates  any 
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significant  difference  between  oceanic  and  continental  attenuation.  To  derive 
the  necessary  regression  equation,  manipulate  equation  (3)  to  get: 


Separate  the  right-hand  side  into  oceanic  and  continental  terms  thus: 


where  r + r = 2irr  , the  total  circumference  of  the  earth.  We  set  T = 20  second 
o c e 

U =3.8  km/sec,  U = 3.0  km/ sec  and  solve  for  the  regression  conef f icients 

A°_1  A _1  C 

Q and  Q by  least  squares  procedure.  The  results  were  inverted  to  get 
6 = 527  and  Q * 341.  A statistical  test  for  identify  of  Q and  Q , assuming 

A0-1  A -1  C -I0  C 

Q - and  Q to  be  random  normal  estimators  of  the  same  Q , showed  that  the 

0 C A -I  A — I A A 

difference  (Q£  - Qq  ),  and  thus  the  difference  between  Qc  and  Qq,  was 

significant  at  only  the  59%  confidence  level.  Use  of  corrected  A^/A^  ratios,  as 

explained  abova,  results  in  Qq  * 618  and  * 406  with  the  difference  being 

significant  at  only  the  60%  confidence  level.  We  cannot  then  reject  the 

possibility  that  Q = Q . Gutenberg  (1945)  and  Tsai  and  Aki  (1969)  report 

o c 

Q < 0 . but  in  these  studies  the  data  is  insufficient  also  to  establish  the 
o c 

relation  between  continental  and  oceanic  attenuation  with  high  confidence. 

The  main  impediment  to  establishing  this  relation  to  a high  degree  of 
confidence  with  our  data  is  that  the  percentages  of  oceanic  versus  continental 
path  as  seen  in  Table  III  are  not  greatly  different,  ranging  from  fifty  to 
eighty  percent  oceanic,  so  that  very  little  of  the  data  scatter  can  be  accounted 
for  by  differences  in  Q between  the  two  structures.  Also,  we  know  that  the 
Rayleigh  waves  examined  will  certainly  travel  paths  different  from  the  great 
circle  path  through  the  epicenter  and  station  and  that  thus  cur  estimated 
percentage  of  oceanic  path  in  Table  III  will  be  a rough  approximation,  possibly 
being  quite  erroneous  in  some  cases. 
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SUMMARY 


•: 


Use  of  large,  carefully  selected  earthquakes,  the  HGLP  network  of  seismic 
stations,  and  simple  digital  processing  of  seismograms  enabled  us  to  routinely 
detect  shorter  periods  of  the  phase.  Amplitude  ratios  of  these  and  asso- 
ciated R^  phases  at  20  seconds'  period  provided  Rayleigh-wave  Q estimates  for 
29  different,  full  great-circle  paths  on  the  globe.  These  empirical  Q 
estimates  include  the  effects  of  reflections,  refractions,  focusing  - 
defocusing,  mode  conversion,  and  scattering  as  well  as  intrinsic  absorption; 
therefore  they  are  relevant  to  determining  standardized  Rayleigh-wave  distance- 

correction  terms  for  calculation  of  M of  worldwide  seismic  events.  Data  of 

s 

this  paper  resulted  in  an  overall  Q of  498,  a value  which  would  predict 
diminution  of  20-second  Rayleigh  waves  to  be  significantly  less  than  that 
implied  by  the  accepted  -1.66  logA  relation  for  amplitude.  Our  data  was 
insufficient  to  show  whether  there  is  a difference  between  oceanic  and 
continental  paths  in  regard  to  attentuation  of  20-second  Rayleigh  waves. 

Although  our  interest  lay  in  the  20-second  periods,  the  data  for  other 
periods  such  as  40  seconds  could  easily  be  gleaned  by  reexamining  the 
processed  seismograms.  In  fact,  the  entire  dispersed  R^  phase  was  vividly 
clear  in  some  cases  after  application  of  the  time-varying  processor  for 
Rayleigh-wave  enhancement;  and  some  features  of  interest  for  further  studies 
appeared,  such  as  multipath  arrivals  for  60-80  second  Airy  phases. 
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